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We introduce a spectral density functional theory which can be used to compute energetics and 
spectra of real strongly-correlated materials using methods, algorithms and computer programs of 
the electronic structure theory of solids. The approach considers the total free energy of a system as 
a functional of a local electronic Green function which is probed in the region of interest. Since we 
have a variety of notions of locality in our formulation, our method is manifestly basis-set dependent. 
However, it produces the exact total energy and local excitational spectrum provided that the exact 
functional is extremized. The self-energy of the theory appears as an auxiliary mass operator 
similar to the introduction of the ground-state Kohn-Sham potential in density functional theory. 
It is automatically short-ranged in the same region of Hilbert space which defines the local Green 
function. We exploit this property to find good approximations to the functional. For example, if 
electronic self-energy is known to be local in some portion of Hilbert space, a good approximation 
to the functional is provided by the corresponding local dynamical mean-field theory. A simplified 
implementation of the theory is described based on the linear muffin-tin orbital method widely used 
in electronic strucure calculations. We demonstrate the power of the approach on the long-standing 
problem of the anomalous volume expansion of metallic plutonium. 

PACS numbers: 71.20.-b, 71.27.+a,75.30.-m 



I. INTRODUCTION 



Strongly correlated electron systems display remarkably interesting and puzzling phenomena, such as high- 
temperature superconductivity, colossal magnetoresistance, heavy fermion behavior, huge volume expansions and 
collapses to name a few. These properties need to be explored with modern theoretical methods. Unfortunately, the 
strongly correlated systems are complex materials with electrons occupying active 3d, 4f or 5f orbitals, (and sometimes 
p orbitals as in many organic compounds and in Bucky-balls-based systems). Here, the excitational spectra over a 
wide range of temperatures and frequencies cannot be described in terms of well-defined quasiparticles. Therefore, the 
design of computational methods and algorithms for quantitative description of strongly correlated materials is a great 
intellectual challenge, and an enormous amount of work has addressed this problem in the pastLSiiiiSiSiLSiSiiSiiiiiS. 

At the heart of the strong-correlation problem is the competition between localization and delocalization, i.e. 
between the kinetic energy and the electron-electron interactions. When the overlap of the electron orbitals among 
themselves is large, a wave-like description of the electron is natural and sufficient. Fermi-liquid theory explains 
why in a wide range of energies systems, such as alkali and noble metals, behave as weakly interacting fermions, 
i.e. they have a Fermi surface, linear specific heat and a constant magnetic susceptibility. The one-electron spectra 
form quasi-particles and quasi-hole bands and the one-electron spectral functions show delta-functions like peaks 
corresponding to the one-electron excitations. We have powerful quantitative techniques such as the density functional 
theory (DFT) in the local density and generalized gradient approximation (LDA and GGA), for computing ground 
state properties^. These techniques can be successfully used as starting points for perturbative computation of one- 
electron spectra, for example using the GW method^. They have also been successfully used to compute the strength 
of the electron-phonon coupling and the resistivity of simple metalsi^. 

When the electrons are very far apart, a real-space description becomes valid. A solid is viewed as a regular array 
of atoms where each element binds an integer number of electrons. These atoms carry spin and orbital quantum 
numbers giving rise to a natural spin and orbital degeneracy. Transport occurs with the creation of vacancies and 
doubly occupied sites. Atomic physics calculations together with perturbation theory around the atomic limit allow 
us to derive accurate spin-orbital Hamiltonians. The one-electron spectrum of the Mott insulators is composed 
of atomic excitations which are broaden to form bands that have no single-particle character. The one-electron 
Green functions show at least two pole-like features known as the Hubbard bandsi^, and the wave functions have an 
atomic-like character, and hence require a many-body description. 

The scientific frontier one would like to explore is a category of materials which falls in between the atomic and 
band limits. These systems require both a real space and a momentum space description. To treat these systems one 
needs a many-body technique which is able to treat Kohn-Sham bands and Hubbard bands on the same footing, and 
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which is able to interpolate between well separated and well overlapping atomic orbitals. The solutions of many-body 
equations have to be carried out on the level of the Green functions which contain necessary information about the 
total energy and the spectrum of the solid. 

The development of such techniques has a long history in condensed matter physics. Studies of strongly corre- 
lated systems have traditionally focused on model Hamiltonians using techniques such as diagrammatic methods^, 
Quantum-Monte Carlo simulations^, exact diagonalizations for finite-size clusters^, density matrix renormalization 
group methods^ and so on. Model Hamiltonians are usually written for a given solid-state system based on physical 
grounds. In the electronic-structure community, the developments of LDA+U^ and self-interaction corrected (SIC)^ 
methods , many-body perturbative approaches based on GW and its extensions^ , as well as time-dependent version 
of the density functional theorji^ have been carried out. Some of these techniques are already much more compli- 
cated and time-consuming comparing to the standard LDA based algorithms, and the real exploration of materials 
is frequently performed by its simplified versions by utilizing such, e.g., approximations as plasmon-pole form for the 
dielectric functioni^, omitting self-consistency within GW^ or assuming locality of the GW self-energjii&. 

In general, diagrammatic methods are most accurate if there is a small parameter in the calculation, say, the ratio 
of the on-site Coulomb interaction U to the band width W. This does not permit the exploration of real strongly 
correlated situations, i.e. when U/W ~ 1. Systems near Mott transition is one of such examples, where strongly 
renormalized quasiparticles and atomic-like excitations exist simultaneously. In these situations, self-consistent meth- 
ods based on the dynamical mean-field based theory (DMFT}i2', and its cluster generalizations such as dynamical 
cluster approximation (DCA)^^, or cellular dynamical mean field theory (C-DMFT)i2*iS, are the minimal many body 
techniques which have to be employed for exploring real materials. 

Thus, a combination of the DMFT based methods with the electronic structure techniques is promising, because a 
realistic material-specific description where the strength of correlation effects is not known a priori can be achieved. 
This work is in its beginning stages of development but seems to have a success. The development was startedSS 
by introducing so-called LDA-I-DMFT method and applying it to the photoemission spectrum of Lai_a:Sra:Ti03. 
Near Mott transition, this system shows a number of features incompatible with the one-electron descriptionSi. The 
LDA-I--I- method^^ has been discussed, and the electronic structure of Fe has been shown to be in better agreement 
with experiment than the one based on LDA. The photoemission spectrum near the Mott transition in V2O3 has 
been studied^^ , as well as issues connected to the finite temperature magnetism of Fe and Ni were explored24. LDA 
-l-DMFT was recently generalized to allow computations of optical properties of strongly correlated materials^^. 
Further combinations of the DMFT and GW methods have been proposed^^'^^'^^ and a simplified implementation to 
Ni has been carried ouli^i. 

Sometimes the LDA-I-DMFT method-'^i omits full self-consistency. In this case the approach consists in deriving 
a model Hamiltonian with parameters such as the hopping integrals and the Coulomb interaction matrix elements 
extracted from an LDA calculation. Tight-binding fits to the LDA energy bands or angular momentum resolved LDA 
densities of states for the electrons which are believed to be correlated are performed. Constrained density functional 
theorjiSS is used to find the screened on-site Coulomb U and exchange parameter J. This information is used in the 
downfolded model Hamiltonian with only active degrees of freedom to explore the consequences of correlations. Such 
technique is useful, since it allows us to study real materials already at the present stage of development. A more 
ambitious goal is to build a general method which treats all bands and all electrons on the same footing, determines 
both hoppings and interactions internally using a fully self-consistent procedure, and accesses both energetics and 
spectra of correlated materials. 

Several ideas to provide a theoretical underpinning to these efforts have been proposed. The eff'ective action ap- 
proach to strongly correlated systems has been used to give realistic DMFT an exact functional formulation^^. Approx- 
imations to the exact functional by performing truncations of the Baym-Kadanoff functional have been discussed^. 
Simultaneous treatment of the density and the local Green function in the functional formulation has been proposedi^. 
Total energy calculations using LDA-I-DMFT have recently appeared in the literaturc'^^''^^''^^'''"'. DMFT corrections 
have been calculated and added to the LDA total energy in order to explain the isostructural volume collapse tran- 
sition in Ce^i. Fully self-consistent calculations of charge density, excitation spectrum and total energy of the S 
phase of metallic Plutonium have been carried out to address the problem of its anomalous volume expansion"^^^. The 
extensions of the method to compute phonon spectra of correlated systems with the applications to Mott insulators^^ 
and high-temperature phases of Ptf^- have been also recently developed. 

In this paper we discuss the details of this unified approach which computes both total energies and spectra of 
materials with strong correlations and present our applications for Pu. We utilize the effective action free energy 
approach to strongly correlated systems22i2Si and write down the functional of the local Green function. Thus, a 
spectral density functional theory (SDFT) is obtained. It can be used to explore strongly correlated materials from 
ab inito grounds provided useful approximations exist to the spectral density functional. One of such approximations 
is described here, which we refer to as a local dynamical mean field approximation. It is based on extended— and 
clusteriiSiiSiiS versions of the dynamical mean-field theory introduced in connection with the model-Hamiltonian 
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approach^''. 

Implementation of the theory can be carried out on the basis of the energy-dependent analog for the one-particle 
wave functions. These are useful for practical calculations in the same way as Kohn-Sham particles are used in density 
functional based calculations. The spectral density functional theory in its local dynamical mean field approximation, 
requires a self-consistent solution of the Dyson equations coupled to the solution of the Anderson impurity model2& 
either on a single siteiS or on a clustenjiiii^. Since it is the most time-consuming part of all DMFT algorithms, we are 
carrying out a simplified implementation of it based on a slave boson Gtuzwiller«Si^Si^ and Hubbard iMiM methods. 
This is described in detail in a separate publication—. We illustrate the applicability of the method addressing the 
problem of (5— Pu. Various aspects of the present work have appeared alread}ii2i2S. 

Our paper is organized as follows. In Section II we describe the spectral density functional theory and discuss local 
dynamical mean field approximation which summarizes the ideas of cluster and extended^ versions of the DMFT. 
We show that such techniques as LDA+DMPT—, LDA+U2., and local GWiSii^ methods are naturally seen within 
the present method. Section III describes our implementation of the theory based on the energy-resolved one-particle 
description^ii and linear-mufhn-tin orbital method^2iiiii for electronic structure calculation. Section IV discusses 
application of the method to the volume expansion of Pu. Section V is the conclusion. 

II. SPECTRAL DENSITY FUNCTIONAL THEORY 

Here we discuss the basic postulates and approximations of spectral density functional theory. The central quantity 
of our formulation is a "local" Green function Gioc{j^,r' , z), i.e. a part of the exact electronic Green function which 
we are interested to compute. This is by itself arbitrary since we can probe the Green function in a portion of a 
certain space such, e.g., as reciprocal space or real space. These are the most transparent forms where the local Green 
function can be defined. We can also probe the Green function in a portion of the Hilbert space. If a function can be 
expanded in some basis set {x^} 

F(r,r',z) = ^Xe(r)%'WxHO (1) 

our interest can, e.g, be associated with diagonal elements of the matrix F^^i{z). 

As we see, the locality is a basis set dependent property. Nevertheless, it is a very useful property because a most 
economical description of the function can be achieved. This is true when the basis set which leads to such description 
of the function is known. The choice of the appropriate Hilbert space is therefore crucial if we would like to find 
an optimal description of the system with the accuracy proportional to the computational cost. In spectral density 
functional theory that has a meaning of finding good approximations to the functional. Therefore we always rely on a 
physical intuition when choosing a particular representation which should be tailored to a specific physical problem. 

At the beginning we formulate spectral density functional theory in completely real space but keep in mind that 
such formulation is not unique. Thus, we are interested in finding a part of the electronic Green function restricted 
within a certain cluster area. Due to translational invariance of the Green function on the original lattice given by 
primitive translations {R}, i.e. G(r + R, r' + R, z) = G(r, r',z), it is always sufficient to consider r lying within a 
primitive unit cell flc positioned at R = 0. Thus, r' travels within some area ilioc centered at R = 0. We set the 
local Green function to be the exact Green function G'(r,r',z) within a given cluster Qioc and zero outside. In other 
words, 

Gioc{r,r',z) = G{r,r',z)eioc{r,r') (2) 

where the theta function is a unity when vector r e f^c, r' G ^loc and zero otherwise. It is schematically illustrated on 
Fig. nj This construction can be translationally continued onto entire lattice by enforcing the property 9ioc (r + R, r' + 

We will now discuss the free energy of a system as a functional of the local Green function. 

A. Functional of Local Green Function 

We consider full many-body Hamiltonian describing the electrons moving in the periodic ionic potential Vext (x) = 
Vext{r)S{T) and interacting among themselves according to the Coulomb law: vc{x — x') — e^/|r — r'\S{T — t') [we 
use imaginary time-frequency formalism, where x = (r,T)]. This is the formal starting point of our all-electron 
first-principles calculation. So, the theory of everything is summarized in the action S:. 
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FIG. 1: Illustration of the area in real space where the local Green function Gioc is defined. Note that r is restricted by the 
unit cell at the origin while r' and r" travel within the crystal. 



S = J dxi;+{x)[dr - V' + Ve^tixMix) 

+ ^ I dxdx'il;'^{x)ip^{x')vc{x ~ x')tl;{x)il;{x') 
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(atomic Rydberg units, h = l,me = 1/2, are used throughout). We will ignore relativistic effects in this action 
for simplicity but considering our applications to Pu, these effects will be included later in the implementation. In 
addition, the effects of electron-phonon interaction will not be considered. 

We will take the effective action functional approach to describe our correlated system^. The approach allows to 
obtain the free energy of the solid from a functional F evaluated at its stationary point. The main question is the 
choice of the variable of the functional which is to be extremized. This question is highly non-trivial because the 
exact form of the functional is unknown and the usefulness of the approach, depends on our ability to construct good 
approximations to it, which in turn depends on the choice of variables. The Baym-Kadanoff (BK) Green function 
theory considers exact Green function G{x,x') = —{Tr'<p{x)'ip^{x')) as a variable, i.e. Tbk[G]. Density functional 
theory considers density p{y) = G(r,r,T = 0) of the solid as a variable, i.e. Tupt[p]- Spectral density functional 
theory will consider local Green function Gioc{x,x') = G{x, x')9ioci'r,r') as a variable, i.e. Tsdf[Gioc]- 

Notice on the variety of choices we can make, in particular in the functional Tsdf[Gioc] since the definition of 
the locality is up to us. The usefulness of a given choice is dictated by the existence of good approximations to 
the functional, as, for example, the usefulness of the DFT is the result of the existence of the LDA or GGA, which 
are excellent approximations for weakly correlated systems. Here we will argue that the usefulness of SDFT is the 
existence of the local dynamical mean field approximations. 

Any of the discussed functionals can be obtained by a Legendre transform of the effective action. The electronic 
Green function of a system can be obtained by probing the system by a source field and monitoring the response. To 
obtain Tbk[G] we probe the system with time-dependent two-variable source field J{x, x') or its imaginary frequency 
transform J(r, r', iuj) defined in all space. If we restrict our consideration to saddle point solutions periodic on the 
original lattice, we can assume that the field obeys the periodicity criteria J(r + R, r' + R, iw) — J{r,r',iuj). This 
restricts the electronic Green function to be invariant under lattice translations. In order to obtain a theory based 
on the density as a physical variable, we probe the system with a static periodical field J{r)S{T). This delivers^^i^Si^ 
the density functional theory Tdft[p]- In order to obtain Tsdf[Gioc\ we will probe the system with a local field 
Jioc{x,x') = Jioc{r,r\T - t') restricted by 6'/oc(r,r'). 

Introduction of the time dependent local source Jioc{x,x') modifies the action of the system as follows 



S'[Jioc] = S + J dxdx Jioc{x,x')Tp{x)Tp^{x) 



(4) 
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Due to translational invariance, the integral over r variable here is the same for any unit-cell flc and the integral 
over r' should be restricted by the area where Jioc ^ 0, i.e. by the cluster area ri;oc -The average of the operator 
^{x)^'^ {x') probes the local Green function which is precisely defined by expression The partition function Z, or 
equivalently the free energy of the system F, becomes a functional of the auxiliary source field 

Z[Jioc] = exp(-F[J,oc]) = / i5[^+V]e^'^'f'''°^' (5) 



The effective action for the local Green function, i.e., spectral density functional, is obtained as the Legendre transform 
of F with respect to the local Green function Gioc{x, x'), i.e. 

^sdf[Gioc] = F[Jioc] — TrJiocGioc (6) 
where we use the compact notation TtJiocGioc for the integrals 

TtJiocGioc= dxdx'Jioc{x,x')Giocix',x)^'Y^ drdr'Jioc{r,r' ,iuj)Giocir' ,r,iu)) (7) 

Using the condition: Jioc = -^ST sdf / SGioc to eliminate Jioc in lO in favor of the local Green function we finally 
obtain the functional of the local Green function alone. 

The source field sets the degree of locality of the object of interest. Considering its definition by expanding the 
cluster till entire solid, we obtain the Baym-Kadanoff functional which determines the Green function in all space. 
Shrinking its definition to a singe point r and assuming its frequency (time) independence, i.e. J(t)5(v — y')5{t — r'), 
we obtain density functional theory. In its extremum, all functionals always reach the total free energy of the system 
regardless the choice of the variable. This situation is similai4&to classical thermodynamics where the thermodynamic 
potential is either the Helmholtz free energy, or the Gibs free energy or the entalpy depending on which variables, 
temperature, pressure, volume are used. Note also that due to assumed time-dependence of the source field, away 
from the extremum the Green function functionals cannot be interpreted as energies. 

Having repeated a formal derivation of the cxistcnce^S of the functional T sY)p\Gioc\ as well as of the functionals 
rB_K-[G'] and Tuft[p\ we now come to the problem of writing separately various contributions to it. This development 
parallels the well known decomposition of the total energy into kinetic energy of a non interacting system, potential 
energy, Hartree energy and exchange-correlation energy. The strategy consists in performing an expansion of the 
functional in powers of the charge of the electron22iiSiiiii2ii£. The lowest order term is the kinetic part of the action, 
and the energy associated with the external potential Vext ■ In the Baym Kadanoff Green function theory this term 
has the form |(3J: 

KBK[G] = Tr\nG-Ti-[GQ^ -G-'^]G (8) 
The Go(r, r', icu) is the non-interacting Green function, which is given by 

Go\r,r',iLj) = 5{r~r')[iLj + fi + V^-Vextir)] (9) 
<^(r-r') = j dv"G^\v,v",zuj)G^{v",v'M (10) 

where is a chemical potential. Note that since finite temperature formulation is adopted we did not obtain simply 
Kbk[G] = Tr( — + Vext)G but also have got all entropy based contributions. 

Let us now turn to the density functional theory. In principle, it does not have a closed formula to describe fully 
interacting kinetic energy as the density functional. However, it solves this problem by introducing a non-interacting 
part of the kinetic energy. It is described by its own Green function Gks{^, r', it^), which is related to the Kohn-Sham 
(KS) representation. An auxiliary set of non-interacting particles is introduced which is used to mimic the density 
of the system. These particles move in some effective one-particle Kohn-Sham potential V^f / (r) — V^xt (r) + Vint (r) . 
This potential is chosen merely to reproduce the density and does not have any other physical meaning at this 
point. The Kohn-Sham Green function is defined in the entire space by the relation G'^^(r, r', iuj) = Gq ^(r, r', iu)) — 
Vi„f(r)(5(r — r'), where Vinti^) is adjusted so that the density of the system p(r) can be found from G'xs(rj r'j «t^)- 
Since the exact Green function G and the local Green function Gioc can be also used to find the density, we can write 
a general relationship: 

p(r) = T ^ GKsir, r, iiu)e'^°+ - T ^ G(r, r, iuj)e'^"+ = T ^ Gioc{r, r, zLu)e'^°+ (11) 
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where the sum over iuj assumes the summation on the Matsubara axis at given temperature T . With the introduction 
of Gks the non-interacting kinetic portion of the action plus the energy related to Vf^^t can be written in complete 
analogy with (jHJ as follows 

Kdft[Gks] = TrlnG/f5 - Tr[Go ^ - G-k\]Gks (12) 
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FIG. 2: Relationship between various Green functions in spectral density functional theory: exact Green function G, local 
Green function Gioc and auxiliary Green function Q are the same in a certain region of space of our interest. They are all 
different outside this area, where the local Green function is zero by definiton. 

In order to describe the different contributions to the thermodynamical potential in the spectral density functional 
theory, we introduce a notion of the energy-dependent analog of Kohn-Sham representation. These auxiliary par- 
ticles are interacting so that they will describe not only the density but also a local part of the Green function of 
the system, and will feel a frequency dependent potential. The latter is a field described by some effective mass 
operator 7V(e//(r, r', ioj) = Vext{i^)S{r — r') + Mint{Tc, r', i^jj)- We now introduce an auxihary Green function Q{r, r', iuj) 
connected to our new "interacting Kohn-Sham" particles so that it is defined in the entire space by the relationship 
Q~^{r,r' ,iuj) = GQ^{r,r' ,iuj) — A1i,it(r, r', zw). Thus, A^i„t(r, r', iw) is a function which has the same range as the 
source that we introduce: it is adjusted until the auxiliary Q(r,r' ,iuj) coincides with the local Green function inside 
the area restricted by 9ioc{r,r'), i.e 

Giocir, r', liu) = g{r, r', H^ioc(r, r') (13) 

We illustrate the relationship between all introduced Green functions in Fig. [3 Note that Q{r,r',iuj) also delivers 
the exact density of the system. With the help of Q the kinetic term in the spectral density functional theory can be 
represented as follows 

KsDF[g]=Trliig -Tr[Go^ -g'']g (14) 

Since Gks is a functional of p, DFT considers the density functional as the functional of Kohn-Sham wave functions, 
i.e. as ToftIGks]- Similarly, since g is a functional of Gioc, it is very useful to view the spectral density functional 
^SDF as a functional of g: 

TsDFig] -^Trlng -TiiG^^ - g-']g + <Psdf[Gioc] (15) 

where the unknown interaction part of the free energy ^sdf[Gioc\ is the functional of Gioc- If the Hartree term is 
explicitly extracted, this functional can be represented as 



<i>SDF[Gloc] = Eh[p] + ^dfIGioc] 



(16) 
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where Eh[p] is the Hartree energy depending only on the density of the system, and where ^s'dfI^Ioc] is the exchange- 
correlation part of the free energy. Notice that the density of the system can be obtained via Gioc or Q, therefore the 
Hartree term can be also viewed as a functional of Gioc or Q. Notice also, that since the kinetic energies (jHl, H12II . 
(|14|l are defined differently in all theories, the interaction energies ^sdf[Gioc\, ^bk[G], ^dft[p] are also different. 
The stationarity of the spectral density functional can be examined with respect to G 

^^'^^ (17) 



6Q(r, r', iuj) 

similar to the stationarity conditions for Tbk[G] and Tdpt[Gks 

STbk 



SG{r,r',iu;) 



= (18) 
= (19) 



6GKsir,r',iu;) 

This leads to the equations for the corresponding Green functions in all theories: 

g-^ (r, r',ioj) = Go\r, r', iuj) - M,nt{r, r', icu) (20) 

as well as 

G~\r,r',iLu) = Go\r,r',iLu)~j:^ntir,r',iu;) (21) 
G^\,(r,r',zc^) = Go\r,r',tLo)~V,nt{r)S{r-r') (22) 

By using Q for Gq ^ and by multiplying both parts by the corresponding Green functions we obtain familiar Dyson 
equations 

[-V^ + Ve^tir) -IUJ- p\Q{v, r', iuj) + ^ dr"M„uir, r",iuj)g{r", r', iuj) = ,5(r - r') (23) 

and 



[-V' + K.t(r)-«w-/i]G(r,r',zc^)+ ^ dr"I],„t(r, r", zc^)G(r", r', icj) = ,5(r - r') (24) 

[-V^ + Ve^t{r)-iu;-ti]GKsir,r',iLu) + Vrnt{r)GKs{r",r',iu;) = d{r-r') (25) 

The stationarity condition brings the definition of the auxiliary mass operator A4mi(r, r', ilu) which is the variational 
derivative of the interaction free energy with respect to the local Green function: 

, , . ^ S^sdf[Gioc] S^sdfIG iOC\ /) / /\ /rtr'\ 
dt;(r',r, iw) SGiocir',r,tu) 

It plays the role of the effective self-energy which is short-ranged (local) in the space. The corresponding expressions 
hold for the interaction parts of the exact self-energy of the electron Si„t(r, r', ioj) and for the interaction part of the 
Kohn-Sham potential T^mt(r). 

E.„.(r,r',z.) = (27) 

V.ndr)Sir - r') ^ J^-^-^^l . = ^^^^.(r - r') (28) 

If the external potential is added to these quantities we obtain total effective self-energies/potentials of the SDF, BK 
and DF theories: Ale//(r, r', zcj), I]e//(r, r', iw), Veff{r) respectively. If the Hartree potential Vff(r) is separated we 
obtain the exchange-correlation parts: A4xciJ^,r' ,iLj), Sa;c(r, r', iw), V^c(r)- 

Note that strictly speaking the substitution of variables, Gks vs. p, in the density functional as well as the 
substitution of variables, Q vs. G/oc, in the spectral density functional is only possible under the assumption of the 
so-called y-representability (or M-representability), i.e. the existence of such effective potential (mass operator) 
which can be used to construct the exact density (local Green function) of the system via the non-interacting Kohn- 
Sham particles of the DFT or its energy-dependent generalization in SDFT. 
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Note also that the effective mass-operator of spectral density functional theory is local by construction, i.e. it is non- 
zero only within the cluster area ^loc restricted by 9ioc(j, r'). It is an auxiliary object which cannot be identified with 
the exact self-energy of the electron Ee//(r, r', icu). This is similar to the observation that the Kohn-Sham potential of 
the DFT cannot be associated with the exact self-energy as well. Nevertheless, the SDFT always delivers local Green 
function and the total free energy exactly (at least in principle) as long as the exact functional is used. In the limit 
when the exact self-energy of the electron is indeed localized within ilioc, the SDFT becomes the Baym-Kadanoff 
functional which delivers the full Green function of the system, i.e. we can immediately identify A4ef f{r,r' jiuj) with 
I]e//(r, r', iw) and the poles of Q{r,r',iuj) with exact poles of G{r,r' ,iuj) where the information about both k and 
energy dependence as well as life time of the quasiparticles is contained. We thus see that, at least formally, increasing 
the size of fijoc in the SDF theory leads to a complete description of the many-body system, the situation quite 
different from the DFT which misses such scaling. 

From a conceptual point of view, the spectral density functional approach constitutes a radical departure from 
the DFT philosophy. The saddle-point equation H23I) is the equation for a continuous distribution of spectral weight 
and the obtained local spectral function Gioc can now be identified with the observable local (roughly speaking, 
k-integrated) one-electron spectrum. This is very different from the Kohn-Sham quasiparticles which are the poles 
of Gks not identifiable rigorously with any one-electron excitations. While the SDFT approach is computationally 
more demanding than DFT, it is formulated in terms of observables and gives more information than DFT. 

On one side, spectral density functional can be viewed as approximation or truncation of the full Baym Kadanoff 
theory where <i>si<-[G] is approximated by ^sdf[Gioc] by restricting G to Gioc^^ and the kinetic functional Kbk [G] 
and Ksdf[G] are thought to be the same. Such restriction will automatically generate a short-ranged self-energy in 
the theory. This is similar to the interpretation of DFT as approximation [G] = ^dft[p], Kbk [G]=Ksjjp[Gks] 
which would generate the DFT potential as the self-energy. However, SDFT can be thought as a separate theory 
whose manifestly local self-energy is an auxiliary operator introduced to reproduce the local part of the Green function 
of the system, exactly like the Kohn-Sham ground state potential is an auxiliary operator introduced to reproduce 
the density of the electrons in DFT. 

Spectral density functional theory contains the exchange-correlation functional ^sdf[Gioc\- An explicit expression 
for it involving a coupling constant X — integration can be obtained in complete analogy with the Harris- Jones 
formula^*' of density functional theory"*^. One considers TsofIQ, A] at an arbitrary interaction A and expresses 

rsDF[g,e']^rsnF[g,0]+ r\x^^^^^^ (29) 

Jo 9X 

Here the first term is simply the kinetic part KsdfIG] sls given by (|14|l which does not depend on A. The second part 
is thus the unknown functional ^sdf[Gioc\- The derivative with respect to the coupling constant in Q is given by the 
average {x)ip~^ {x')ip{x)^{x')) = Ux{x,x' ,iu;)+ {ip'^ {x)^{x)){ip~^ {x')'il){x')) where Il\{x,x') is the density-density 
correlation function at a given interaction strength A computed in the presence of a source which is A dependent and 
chosen so that the local Greens function of the system is Q. Since {^p~^ {x)ip {x)) = p{r)S{T), we can obtain : 

2 

<i>5DF [G,oc] = Eh [p] r dA^^^^^^ii^ (30) 



Establishing the diagrammatic rules for the functional ^sdf[Gioc] while possibleSS, is not as simple as for the 
functional <&_Bi<-[G]. The latter is formally represented as a sum of two-particle diagrams constructed with G and vc- 
It is well known that instead of expanding <I>bx[G] in powers of the bare interaction vc and G, the functional form can 
be obtained by introducing the dynamically screened Coulomb interaction W{r,r',ioj) as a variable. In the effective 
action formalism^ this was done by introducing an auxiliary Bose variable coupled to the density, which transforms 
the original problem into a problem of electrons interacting with the Bose field. W is the connected correlation 
function of the Bose field. 

Our effective action is now a functional of G, W and of the expectation value of the Bose field. Since the latter 
couples linearly to the density it can be eliminated exactly, a step which generates the Hartree term. After this 
elimination, the functional takes the form 

rBK[G,W] =Ti\nG-Tr[G^^ -G-^]G + <^bk[G,W] (31) 



$bk[G, VF] = Eh[p] - -TyIuW + -Tr[vc^ - W-^]W + ^bk[G,W] (32) 
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The entire theory is viewed as the functional of both G and W. Here, ^ bk[G,W] is the sum of all two-particle 
diagrams constructed with G and W with the exclusion of the Hartree term, which is evaluated with the bare 
Coulomb interaction. An additional stationarity condition STbk/SW — leads to the equation for the screened 
Coulomb interaction W 

W-\r,r',iuj)^v^\r~r')~U{r,r',ioj) (33) 

where the function n(r,r',zcj) = —2d'i>BK/SW{r,r',iuj) is the exact interacting susceptibility of the system, which is 
already discussed in connection with representation (|30|) . 

A similar theory is developed for the local quantities"^*^, and this generalization represents the ideas of extended 
dynamical mean field theory^^^, now viewed as an exact theory. Namely, one constructs an exact functional of the 
local Greens function and the local correlator of the Bose field coupled to the density which can be identified with 
the local part of the dynamically screened interaction. The real-space definition of it is the following 

Wiocir, r', iu) = W{v, r', ic^)0/o,(r, r') (34) 

which is non-zero within a given cluster f2;oc- Note that formally this cluster can be different from the one considered 
to define the local Green function |(2Jl but we will not distinguish between them for simplicity. An auxiliary interaction 
W(r, r', iut) is introduced which is the same as the local part of the exact interaction within non-zero area of dioc{^, r') 

Wiocir, r', tiu) = W(r, r', iu;)9iocir, r') (35) 

The interaction part of the spectral density functional is represented in the form similar to H32|l 

^sdf[Gioc,Wioc] = Eh[p] - ^Ti-lnW +^Tt[vc^ -W-^]W + ^sdf[Gioc,Wioc] (36) 

and the spectral density functional is viewed as a functional Tsdf[Gioc,Wioc\ or alternatively as a functional 
^SDf[G, VV]. sdf[Gioc, Wioc] is formally not a sum of two-particle diagrams constructed with Gioc and Wioc, but 
in principle a more complicated diagrammatic expression can be derived. Alternatively, a more explicit expression 
involving a coupling constant integration can be given. Examining stationarity 5T sdf / 5W = yields a saddle-point 
equation for W(r, r', iui) 

yy"^(r,r',iw) ^v^^iv -v') - V{v,v' ,iuj) (37) 
where the effective susceptibility of the system is the variational derivative 



P(r, V ,iuS) = — — — = , -^6';oc(r, r ) 38 

d>V(r',r,ia;) bWioc{'!^' ,v,iuS) 

Notice again a set of parallel observations for V as for Meff, Eq- (|26|l . The effective susceptibility of spectral density 
functional theory is local by construction, i.e. it is non-zero only within the cluster restricted by ^ioc(r,r')- Formally, 
it is an auxiliary object and cannot be identified with the exact susceptibility of the electronic system Il{r , r' , iuj) . 
However, if the exact susceptibility H(r,r',iaj) is sufficiently localized, this identification becomes possible. If cluster 
rijoc includes physical area of localization, we can immediately identify V{r,r', iuj) with Il{r,r',iuj) and W{r,r',iuj) 
with W{r,r',iuj) in all space. However, both W and W are always the same within ftioc regardless its size, as it is 
seen from (|34|) and H35|l . 

At the stationarity point, TsdfIG, W] is the free energy F of the system. If one inserts (|20|) into (|14l) and (|37|) into 
(|36|l we obtain the formula: 

F = Trln^ - TiMeffG + TvVe,tG+EH - ^Trln W + ^TrPW + ^sdf (39) 
Similar formulae hold for the Baym-Kadanoff and density functional theories 

F = Tr In G - TrEe/yG + Tr^^tG + Eh - ^Tr InW + ^TtHW + -^bk (40) 

F = TvIuGks -TrVeffGxs + TrVe^tGKS + '^DFT (41) 

where the first two terms in all expressions (|39|l . H4l)|) . H41|) are interpreted as corresponding kinetic energies, the third 
term is the energy related to the external potential Vext which is in fact TrVextP in all cases. The other terms represent 
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the interaction parts of the free energy. Note that all entropy originated contributions are included in both kinetic 
and interaction parts. If temperature goes to zero, the entropy part disappears and the total energy formulae will be 
recovered. For example, in spectral density functional theory we obtain: 

E = -TrV'e + TrKxtP + Eh + (42) 

We will also discuss this limit later in more details in Section III. 

The SDFT approach is so far not very useful since a tractable expression for the functional form of '^sdf[Gioc\ or 
^ sdf[Gioc, VVioc] has not been given yet. This is quite similar to the unknown exchange-correlation functional of the 
DPT. As we have learned from the developments of the dynamical mean-field methods, a very useful approximation 
exists to access these functionals. This is based on a full many-body solution of a finite-size cluster problem treated 
as an impurity embedded into a bath subjected to a self-consistency condition. Such local dynamical mean field 
theory will be discussed below. 



B. Local Dynamical Mean Field Approximation 

The spectral density functional theory, where an exact functional of certain local quantities is constructed in the 
spirit of Ref. uses effective self-energies and susceptibilities which are local by construction. This property can 
be exploited to find good approximations to the interaction energy functional. For example, if it is a priori known 
that the real electronic self-energy is local in a certain portion of the Hilbert space, a good approximation is the 
corresponding local dynamical mean field theory obtained for example by a restriction or truncation of the full Baym- 
Kadanoff functional or its generalization to use W and G as natural variables, to local quantities in the spirit of Ref. 

M 

The local DMFT approximates the functional ^sdf[Gioc] (or sdf[Gioc,Wioc]) by the sum of all two-particle 
diagrams evaluated with Gioc and the bare Coulomb interaction vc (or screened local interaction Wioc)- In other 
words, the functional dependence of the interaction part ^^^[G] in the Baym-Kadanoff functional for which the 
diagrammatic rules exist is now restricted by Gioc and is used as an approximation to ^sdf[Gioc\, i-e. ^sdf[Gioc\ — 
$ BK [Gioc] ■ Obviously that the variational derivative of such restricted functional will generate the self-energy confined 
in the same area as the local Green function itself. 

Remarkably the summation over all local diagrams can be performed exactly via introduction of an auxiliary 
quantum impurity model subjected to a self-consistency conditioniSi^. If this impurity is considered as a cluster C, 
the cellular DMFT (C-DMFT) can be used which breaks the translational invariance of the lattice to obtain accurate 
estimates of the self energies. The C-DMFT approximation, can also be motivated using the cavity construction. 
The solid should be separated onto large cells which circumscribe the areas ri/oc- Considering the effective action S, 
Eq. Q, the integration volume is separated onto the cellular area flc and the rest bath area — = ^bath- The 
action is now represented as the action of the cluster cell, Qc plus the action of the bath, flbath, plus the interaction 
between those two. We are interested in the local effective action Sc of the cluster degrees of freedom only, which is 
obtained conceptually by integrating out the bath in the functional integral: 

^exph5c] = ^/ drdr'exp[-5] (43) 

where Zc and Z are the corresponding partition functions. This integration is carried out approximately, keeping 
only a charge-charge interaction as quartic terms and neglecting all the higher order terms generated in this process 
to arrive to a cavity action of the formiSiSSiSSiS^: 

Sc = - j dx4i^{x)QQ^{x,x)i}j{x) 

+ ^ J dxdx'iP+{x)i:+{x')Voix,x')^{x)ip{x') (44) 

where the integration over the spatial variables is performed over Qc- Here Goi^, x') or its Fourier transform Qo{y, r', iuo) 
is identified as the bath Green function appeared in the Dyson equation for the local mass operator A^i„t(r, r', zcj) 
and for the local Green function Gioc{v,v' ^iui) of the cluster, and Vo(r,r',za;) is the "bath interaction" appeared in 
the Dyson equation for the local susceptibility ■p(r,r',?a;) and local interaction Wioc(r, r', i.e 

g^\ry,iu) = GY^\{vy,iLo)+MM[r,v',iLo) (45) 

V^\Y,Y\iuj) = W i-^l [y y ,Ujj) +V{Yy .,10J) (46) 
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Note that neither nor Vq can be associated with non-interacting Go and bare interaction vc, respectively. Note 
also that both r and r' indexes in t/o(r, r', and in Vo(r, r', iuj) vary within the cellular area ^Ic- The same should 
be assumed for the local quantities Gzoc(r, r', ioj) and Wioc(r, r', Since these functions are translationally invariant 
on the original lattice, this property can be used to set up these functions within Qc- 

An interesting observation can be made on the role of the impurity model which in the present context appeared 
as an approximate way to extract the self-energy of the lattice using input bath Green function and bath interaction. 
Alternatively, the impurity problem can be thought itself as the model which delivers exact mass operator of the 
spectral density functional^. If the latter is known, there should exist such bath Green function and such bath 
interaction which can be used to reproduce it. In this respect, the local interaction Wioc appeared in our formulation 
can be thought as an exact way to define the local Coulomb repulsion "[/", i.e. such interaction which delivers exact 
local self-energy. 

To summarize, the effective action for the cluster cell (|44|) . the Dyson equations H45|) . (|46l) connecting local and 
bath quantities as well as the original Dyson equations (|20|l . (|37|l constitute a self-consistent set of equations as the 
saddle-point conditions extremizing the spectral density functional Tsdf{G, VV). They combine cellular and extended 
versions of DMFT and represent our philosophy in the ab initio simulation of a strongly correlated system. Since 
Mint and V are unknown at the beginning, the solution of these equations assumes self-consistency. First, assuming 
some initial 7Wi„t,and V the original Dyson equations H20|l . (|37() are used to find Green function Q and screened 
interaction W. Second, the Dyson equations for the local quantities H45|l . (|46|) are used to find Qq, Vq. Third, quantum 
impurity model with the cluster action Sioc after (|44|l is solved by available many-body technique to give new local 
Mint and V. The process is repeated till self-consistency is reached. This is schematically illustrated in Fig. El Note 
here that while single-site impurity problem has a well-defined algorithm to extract the lattice self-energy, this is not 
generally true for the cluster impurity modelaiS. The latter provides the self-energy of the cluster, and an additional 
prescription such as implemented within cellular DMFT or using DCA should be given to construct the self-energy 
of the lattice. 
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FIG. 3: Illustration of self-consistent cycle in spectral density functional theory with local dynamical mean-field approximation: 
both local Green function Gioc and local Coulomb interaction Wioc are iterated. The auxiliary quantities Q and W are used to 
simplify the construction of the functional. 



Unfortunately, writing down the precise functional form for ^sdfY^Ioc^ ^ioc\ or ^ SDF^loc^ Wioc] is still a problem 
because the evaluation of the entropy requires the evaluation of the energy as a function of temperature and an 
additional integration over it. In general, the free energy F ^ E — TS, where E is the total energy and S is the 
entropy. Since Tsup[Q] = KsDp[G] + ^sdf[Gioc\j both energy and entropy terms exist in the kinetic and interaction 
functional. The energy part of Ksdf[Q] = Tr( — -f Vext)Q and the energy part of ^sdf[Gioc,Wioc\ can be 
written explicitly as ^TiMintGioc- The entropy correction is a more difficult one. In principle, it can be evaluated 
by performing calculations of the total energy Esdf[G] = Tr( — + Vext)G + ^^TMintGioc at several temperatures 
and then taking the integraliii: 

SiT) = S{oo)-j^dr^^-^ (47) 

The infinite temperature limit 5(00) for a well defined model Hamiltonian can be worked out. This program was 
implemented for the Hubbard model^^ and for Ga^i. 
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Two well separate problems are now seen. For a given material using the formulae H20|l . H37() . 101, H34() . H45() . H46() 
Cy, W, Gioc, W^/oc7 ^/Oi ^0 should be computed using the methods and algorithms of the electronic structure theory. This 
procedure will in part be described in Section III. With given input Qq and Vq, the solution of the impurity model 
constitutes a well separated problem which can be carried out either using the QMC method or other impurity solver. 
Some of the techniques are discussed in Refs. .ICLllt In Section IV, while applying a simplified version of the theory 
to plutonium, we will briefly describe an impurity solver used in that calculation. A full description of this method 
will appear elsewhere^. 

The described algorithm is quite general, totally ab initio and allows to determine all quantities, such as the one- 
electron local Green functions G/oc and the dynamically screened local interactions Wioc- Unfortunately, its full im- 
plementation is a very challenging project which so far has only been carried out at the level of model HamiltoniansS^. 
There are several simplifications which can be made, however. The screened Coulomb interaction W(r,r',iuj) can 
be treated on different levels of approximations. In many cases used in practical calculations with the LDA+DMFT 
method, this interaction W is assumed to be static and parametrized by a set of some optimally screened on-site 
parameters, such as Hubbard U and exchange J. These parameters can be fixed by constrained density functional 
calculations, extracted from atomic spectra data or adjusted to fit the experiment. Since the described theory can 
perform a search in a constrained space with fixed interaction W, this justifies the use of U and J as input num- 
bers. A more refined approximation, can use a method such as GW to generate an energy-dependent W-' which 
is then treated using extended DMFT^^. Alternatively we can envision that W is already so short ranged that we 
can ignore the EDMFT self consistency condition, and we treat W as Wfi^ix, x') .This leads to performing a partial 
self-consistency with respect to the Green function only. The procedure is reduced to solving Dyson equations 1(20)), 
(|45|l as well as to finding Mint via the solution of the impurity problem. A full self-consistency can finally be restored 
by including a second loop to relax W. 

A methodological comment should be made in order to make contact with the literature of cluster extensions of 
single site DMFT within model Hamiltonians. We adopted a less restrictive notion of locality by defining an effective 
action of the one particle Green function (and of the interaction) whose arguments are in nearby unit cells. This 
maintains the full translation invariance of the lattice. At the level of the exact effective action , this is an exact 
construction, and its extremization will lead to portions of the exact Greens function which obeys causality. Notice 
however that it has been proved recentljiiS that generating approximations to the exact functional by restricting the 
Baym Kadanoff functional to non local Greens functions leads to violations of causality. For this reason, we propose 
to use techniques such as CDMFT which are manifestly causal for the purpose of realizing approximations to the 
local Greens functions. 

Our final general comment concerns the optimal choice of local representation or, precisely, the definition of the 
local Green function. This is because the local dynamical mean-field approximation is likely to be accurate only if 
we know in which portion of the Hilbert space the real electronic self-energy is well localized. Unfortunately, this 
is not known a priori, and in principle, only a full cluster DMFT calculation is capable to provide us some hints in 
attempts to answer this question. However, considerable empirical evidence can be used as a guide for choosing a 
basis for DMFT calculations, and we discuss these issues in the following sections. 

C. Choice of Local Representation 

We have already pointed out that spectral density functional theory is a basis set dependent theory since it probes 
the Green function locally in a certain region determined by a choice of basis functions in the Hilbert space. Provided 
the calculation is exact, the free energy of the system and the local spectral density in that Hilbert space will be 
recovered regardless the choice of it. We have developed the theory assuming that the basis in the Hilbert space is 
indeed the real space which gives us the choice Q for the local Green function, i.e. the part of the real Green function 
restricted by 0/oc(r, r'). While this is most natural choice for the purpose of formulating locality in r and r' variables, 
it is also very useful to discuss a more general choice, connected to some space of orbitals X^(r) which can be used 
to represent all the relevant quantities in our calculation. As we have in mind to utilize sophisticated basis sets of 
modern electronic structure calculations, we will sometimes waive the orthogonality condition and will introduce the 
overlap matrix O^^' — (x^lxf) especially in cases when we discuss a practical implementation of the method. 

We note that the space Xj(r) can in principle be interpreted as the reciprocal space plane wave representation 
Xj(r) = 6*^'^+*^''",^ = k + G with k being the Brillouin zone vector and G being the reciprocal lattice vector. Thus 
the Green function can be probed in the region of the reciprocal space. It can be interpreted as the real space 
representation if x^(r) = — r) where the sums over ^ are interpreted as the integrals over the volume, and 
the locality in this basis is precisely exploited in ((^J. A tremendous transparency of the theory will also arrive if 
we interpret the orbital space {x^} as a general non-orthogonal tight-binding basis set when index ^ combines the 
angular momentum index Im, and the unit cell index R, i.e. x^i^) = Ximi^ ~ R-) = Xa ('^ ~ I^) ■ N^ote that we can add 
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additional degrees of freedom to the index a such, for example, as multiple kappa basis sets of the linear muffin-tin 
orbital based methods, Gaussian decay constants in the Gaussian orbital based methods, and so on. If more than one 
atom per unit cell is considered, index a should be supplemented by the atomic basis position within the unit cell, 
which is currently omitted for simplicity. For spin unrestricted calculations a accumulates the spin index a and the 
orbital space is extended to account for the eigenvectors of the Pauli matrix. 

Let us now introduce the representation for the exact Green function in the localized orbital representation 

G(r, v',zcu)^Y.T. X^(r)G„;3(k, ^io)xf{r') = ^ ^ Xa(r - R)G„^(R - R', iio)x*pir' - R') (48) 

a/3 k a/3 RR' 

Assuming the single-site impurity case, we can separate local and non-local parts Giod^, r', iuj) + Gnon-ioc{'^i r', iw) 
as follows 

G,oc(r, r', zc^) = ^ Gio,^^p{iu:) ^ Xa(r " ^)x}{r' - R) - E Gloc.,o.p{^u:) ^ Xa(r)x^*(r') (49) 

a/3 a/3 k 

where we denoted the site-diagonal matrix elements 5 rw G ap^R. — R^iw) as Giocapii^)- Note that this definition 
is different from the real-space definition (0). For example, |(5J contains the information about the density of the 
system. The formula H49|l does not describe the density since R 7^ R' elements of the matrix Gq/3(R — R',ia;) are 
thrown away. The locality of H49|l is controlled exclusively by the decay of the orbitals Xai'^) ^ ^ function of r, not 
by 0/oc(r,r') 

The local part of the Green function, G;oc(r, r', i^^), which is just defined with respect to the Hilbert space {Xa) 
can be found by developing the corresponding spectral density functional theory. Since the basis set is assumed to be 
fixed, the matrix elements Gioc,ap{i^) appear only as variables of the functional. As above, we introduce an auxiliary 
Green function t/a/3(k, iuj) to deal with kinetic energy counterpart. Stationarity yields the matrix equation: 

G"^^ (k, iu) = g-^ (k, iw) + Mint,af3 {ioj) (50) 

where the non-interacting Green function is the matrix of non-interacting one-electron Hamiltonian 

Go~i^(k, tLo) - (x^ + fi + V'- V;,t Ix^) (51) 

The self-energy Mmt.a/sii'^) is the derivative 6^sDF[Gioc.a[3ii'-^)]/SGioc.ai3{i^) and takes automatically the k- 
independent form. 

While formally exact, this theory would have at least one undesired feature since, for example, the density of the 
system can no longer be found from the definition l|49(l of G;oc(r, r', zcj). As a result the Hartree energy cannot be 
simply recovered. If treated exactly ^sopiGiocafsii^)] should contain the Hartree part. However, we see that the 
theory delivers k-independent Mint.at^ii'-^) including the Hartree term. There seems to be a paradox since modern 
electronic structure methods calculate the matrix element of the Hartree potential within a given basis exactly, i.e. 
(Xa|V//|x^)- The k-dependence is trivial here and is connected to the known k-dependence of the basis functions 
used in the calculation. Therefore, while formulating the spectral density functional theory for electronic structure 
calculation, we need to keep in mind that in many cases, the k-dependence is factorizable and can be brought into 
the theory without a problem. This warns us that the choice of the local Green function has to be done with care so 
that useful approximations to the functional can be worked out. It also shows that in many cases the k-dependence 
is encoded into the orbitals. It is not that non-trivial k-dependence of the self-energy operator, which is connected 
to the fact that Mint (r,r',*w) may be long-range, i.e. decay slowly when r departs from r'. It may very well be 
proportional to 6{r — r') like the LDA potential and still deliver the k-dependence. 

It turn out that the desired k dependence with the choice of the Green function after (|49|l can be quickly reinstated 
if we add the density of the system as another variable to the functional. This is clear since the density is a particular 
case of the local Green function in 10 ) taken at r = r' and summed over iuj. Therefore combination of definition 
(|49|l and p is another, third possibility of defining Gioc- This will allow treatment of all local Hartree-like potentials 
without a problem. Moreover, as we discuss below, this may allow to design better approximations to the functional 
since the Hilbert space treatment of locality is more powerful: it may allow us to treat more long-ranged self-energies 
than the ones restricted by 0;oc(r,r'), and the basis sets can be optimally adjusted to specific self-energies exactly as 
the basis sets used in electronic structure calculations are tailored to the LDA potential. 

We have noted earlier that the mass operator A4i„((r, r', iuj) is an auxiliary object of the spectral density functional 
theory. It has the same meaning as the DFT Kohn-Sham potential: it is local self-mass operator that needs to be 
added to the non-interacting Green function in order to reproduce the local Green function of the system, as the DFT 
potential is added to the non-interacting Green function to reproduce the density of the system. Roughly speaking. 
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SDFT provides the exact energy and exact one-electron density of states which is advantageous compared to the DFT 
which provides the energy and the density only. However, we obtain the full k-dependent one-particle spectra as 
the poles of auxiliary Green function 0{r, r' ,z). Can these poles be interpreted as the exact k-dependent one-electron 
excitations? This question is similar to the question of the DFT: can the Kohn-Sham spectra be interpreted as the 
physical one-electron excitations? To answer both questions we need to know something about exact self-energy of 
the electron. If it is energy-independent, totally local, i.e. proportional to 6{r — r') and well-approximated by the 
DFT potential, the Kohn-Sham spectra represent real one-electron excitations. The exact SDFT waives most of 
the restrictions: if the real self-energy is localized within the area Rioc, the exact SDFT calculation with the cluster 
flioc including Rioc will find the exact k-dependent spectrum. If we pick flioc larger than i?/oc, the SDFT equations 
themselves will choose physical localization area for the self-energy during our self-consistent calculation. However, 
these statements become approximate if we utilize the local dynamical mean field approximation instead of extremizing 
the exact functional. Even if the real self-energy of the electron is sufRciently short-ranged, this approximation will 
introduce some error in the calculation, the situation similar to LDA within DFT. However, the local dynamical mean 
field theory does not necessarily have to be formulated in real space. The assumption of localization for self-energy 
can be done in some portion of the Hilbert space. In that portion of the Hilbert space the cluster impurity model 
needs to be solved. 

The choice of the appropriate Hilbert space, such, e.g., as atomic-like tight-binding basis set is crucial, to obtain 
an economical solution of the impurity model. Let us for simplicity discuss the problem of optimal basis in some 
orthogonal tight-binding (Wannier-like) representation for the electronic self-energy 

I](r, r', zc.) = ^ ^ Xo(r)Sa/3(k, zc.)x^*(r') = E E ~ R)S„/3(R - R', ^u:)x}ir' - R') (52) 

Q/3 k Q/3 BR' 

We can separate our orbital space {Xa} onto the subsets describing hght {xa} ^^'^ heavy {Xa} electrons. Assum- 
ing either off-diagonal terms between them are small or we work with exact Wannier functions, the self-energy 
S(r,r',ia;) can be separated onto contributions from the light, I]i(r, r', zcj), and from the heavy, E// (r, r', iw), elec- 
trons. T,af3{'k,iuj) is expected to be k-dependent but largely ui independent for the hght block, i.e Si(r, r', jw) = 
J2 AB J2k x!4(r)5]AB(k)xB*(r')- The k-dependency here should be well-described by LDA-like approximations, there- 
fore we expect Ei(r,r',iw) ~ T4//(r)(5(r — r'). A different situation is expected for the heavy block where we would 
rely on the result 

E^(r,r',zw) ~ Veff{r)6{r - r') -t- ^ XaWS^.Mx^ (r') (53) 

ab 

The first term here gives the k-dependence coming from an LDA-like potential. It describes the dispersion in the 
heavy band. The second term is the energy dependent correction where site-diagonal approximation R = R' is 
imposed. What is the best choice of the basis to use in connection with '^'^(^{iuj) in l|53() ? Here the decay of the 
orbitals Xai^) ^ function of r is now entirely in charge of the self-energy range. In light of the spectral density 
functional theory, the answer is the following: the local dynamical mean field approximation would work best for such 
basis Xa (i") whose range approximately corresponds to a self-energy localization Rioc of the real electron. Even though 
Rioc is not known a priori, something can be learned about its value based on a substantial empirical evidence. It is, 
for example, known that LDA energy bands when comparing to experiments at first place miss the energy dependent 
T,'^i^{iu}) like corrections. This is the case for bandwidths in transition metals (and also in simple metals), the energy 
gaps of semiconductors, etc. It is also known that many-body based theories work best for massively downfolded 
model Hamiltonians where only active low-energy degrees of freedom at the region around the Fermi level Ep remain. 
The many-body Hamiltonian 

^ = E E ^^%R' + '^•^■l + E E V^77''''''c+^c+^,csR>,,c,n'' (54) 

Q/3 RR' al3iS RR'K'R'" 

with V^p^f^"^'" = J drdr'x*7?(r)x^fl/(r')wc(r - i"')X7fl" (i")Xa-_r"' (i"') assumes the one-electron Hamiltonian hfj^^j^ 
is obtained as a fit to the bands near Ep. This can always be done by long ranged Wannier functions. It is also 
clear that the correlation effects are important at first place for the partially occupied bands since only these bring 
various configurational interactions in the many-body electronic wave functions. For example, the well-known one- 
band Hamiltonian for Cu02 plane of high-Tc materials considers an antibonding combination of Cu2;2_j,2 and O^r^j^ 
orbitals which crosses Ep. Also, the calculations based on the LDA+DMFT method usually obtain reliable results 
when treating only the bands crossing the Fermi level as the correlated one-electron states. This is, for example. 



15 



the case of Pu or our— and previous^ calculation for LaTiOs where t2g three band Hamiltonian is considered. All 
this implies that the range for T,'^i^{iuj) term in (|53|1 should correspond to the properly constructed Wannier orbitals, 
which is fairly long-ranged. What happen if we instead utilize mostly localized representation which, for example, 
can be achieved by tight-binding fits to the energy bands at higher energy scale. For the case of Cu02 this would 
correspond to a three band Hamiltonian with Cu3.2_y2 and Ox.y orbitals treated separately. For LaTiOs system this 
is a Hamiltonian derived from Tit2g and Op orbitals. The answer here can be given as a practical matter of most 
economic way to solve the impurity problem: provided Cu and O levels are well separated, provided both approaches 
use properly downfolded for each case Coulomb interaction matrix elements Va 'p^s^'^ ' provided correlations are 
treated on all orbitals, the final answer should be similar regardless the choice of the basis. A faster algorithm will be 
obtained by treating the one-band Hamiltonian with antibonding Cu2;2_j^2- Ox^y orbital. If indeed the self-energy is 
localized on the scale of the distance between Cu and O, it is clear where the inefficiency of the three-band model 
appears: the second term in (|53|l needs to be extended within the cluster involving both Cu and nearest O sites and 
should involve both Cu and O centered orbitals simply to reach the cluster boundary. In the one-band case this is 
encoded into the decay of the properly constructed Wannier state. 

The previous discussion is merely a conjecture. It does not imply that the localization range for the real self-energy 
of correlated electron at given frequency uo is directly proportional to the size of Wannier states located in the vicinity 
of w + /i. It may very well be that in many cases this range is restricted by a single atom only (atomic sphere of Cu 
in the example above). Clearly more experience can be gained by studying a correlation between the decay of the 
Coulomb matrix element V^^^ ■'^ as a function of R — R' and the obtained matrix E(R — R'jt^^) using a suitable 
cluster DMFT technique. These works are currently being performed and will be reported elsewhere^. The given 
discussion however warns that in general the best choice of the basis for single-site dynamical mean field treatment 
may not be the case of mostly localized representation. In this respect the area restricted by 0/oc(r, r') which is used 
to formulate SDFT in the real space may need to be extended up to a cluster. However, alternative formulation 
with the choice of local Green function after H49|l may be more economical since a single-site approximation may 
still deliver good results. As we have argued, such spectral density functional theory will also need the density of 
the system to complete the definition of local Green function. The local dynamical mean field approximation can be 
applied to the interaction functional ^sdf which is viewed as $sdf[p, G/oc]- This idea is used by the LDA+DMFT 
method described below. 



D. LDA+DMFT Method 

Various methods such as LDA+U^, LDA+DMFT^^ and local GWi2ii& which appeared recently for realistic calcu- 
lations of properties of strongly correlated materials can be naturally understood within spectral density functional 
theory Let us, for example, explore the idea of expressing the energy as the density functional. Local density ap- 
proximation prompts us that a large portion of the exchange-correlation part 'I'a^cip] can be found easily. Indeed, 
the charge density is known to be accurately obtained by the LDA. Why not think of LDA as the most primitive 
impurity solver, which generates manifestly local self-energy with localization radius collapsed to a single r point? It 
is tempting to represent '^sdf[Gioc] = Eh[p] + E^l^^[p\ + ^[Gioc] - '^dc[Gioc], where the new functional ^sdf[Gioc] 
needs in fact to take care of those electrons which are strongly correlated and heavy, thus badly described by LDA. 
Conceptually, that means that the solution of the cluster impurity model for the light electrons is approximated by 
LDA and does not need a frequency resolution for their self-energies. 

Unfortunately, the LDA has no diagrammatic representation, and it is difficult to separate the contributions from 
the light and heavy electrons. The E^^^lp] is a non-linear functional and it already includes the contribution to 
the energy from all orbitals in some average form. Therefore we need to take care of a non-trivial double counting, 
encoded in the functional ^dc[Gio(^- The precise form of the double counting is related to the approximation imposed 
for ^[Gioc\- We postpone this discussion until establishing the connection to the LDA+U method in the following 
subsection. 

The LDA-I-DMFT approximation considers both the density and the local Green function Giocapii^) defined in 
(|49|l as the parameters of the spectral density functional^^ . A further approximation is made to accelerate the solution 
of a single-site impurity model: the functional dependence comes from the subblock of the correlated electrons only. 
If localized orbital representation {Xa\ is utilized, a subspace of the heavy electrons {Xal can be identified. Thus, the 
approximation can be written as ^SDF[Gioc,ab{i'^)] , where Gioc,ab{i^) is the heavy block of the local Green function. 
The double counting correction depends only on the average density of the heavy electrons. Its precise form will 
be discussed below, but for now we assume that (^dc[Gioc] = 4'z5c[^c] with fie = Tj^i^j J2a Gioc.aa{'it^)e^'^^ , where 
index a runs within a correlated Ic shell only. We can write the LDA+DFMT approximation for the interaction energy 
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as follows: 

<^lda+dmft[p, Gioc] = Eh[p] + E^^%] + ^Gioc,ab{iLo)] - ^Dc[nc] (55) 

The kinetic energy part is treated as usual with introducing the auxiliary Green function Q{v, r', ilo). 

The full functional T l£,a+dmft[Q] is considered as a functional of the matrix Qap^H — R',ia;) or its Fourier 
transformed analog ^Q,^(k, iw). The stationarity is examined with respect to ^Q,^(k, iw) and produces the saddle- 
point equation similar to (j2()|l . It has the following matrix form 

^0^1/3 C^'*^) = ^Q/3(k, it^) +7Wmt,Q/3(k, icj) (56) 

where the non-interacting Green function ((SJ is the matrix of non-interacting one-electron Hamiltonian 

G-i^(k, iLo) = {xlV^ + A* + - V,^t\x),) (57) 

The self-energy ^Aint,ap{^,il^) is the variational derivative of '^lda+dmft[p,Gioc]- Its precise form depends on the 
basis set used in the LDA+DMFT calculation. 

In general, it can be split onto several contributions including Hartree, LDA exchange-correlation, DMFT and 
the double-counting correction. In orthogonal tight-binding, both DMFT, M.ab{iw), and double counting, 5abV^'" , 
matrices do not depend on k. These matrices are non-zero within the heavy block only. The Dyson equation H56|) 
can be rewritten by separating from Mint,a0{}^, i^^) the total LDA potential Vlda{y) — Vext{'r) + Vh{y) + V'^^-°^(r): 

The Green function Q^pi)^, ii^) obtained from is used to find Gioc.afiii^) — X^k Qap{}^j i^) which is then used in 
another Dyson equation to compute the bath Green function: 

Q^lbi^^) = GtX^,(,iu) + Mabiiu) (59) 

In Section III we will also describe an accurate procedure to solve the real space form (j45|l of the Dyson equation using 
the LMTO basis set. The LDA+DMFT bath Green function ^o,ab(*'^) is the only essential input to the auxiliary 
impurity model. Thus, the procedure of self-consistency within LDA+DMFT is reduced to the following steps. First, 
some self-energy matrix of the heavy orbitals M.ab{i'^) is guessed. Then, the Dyson equation H56|) is solved in the 
entire Hilbert space and delivers the Green function ^^^(k, ioj). After that, the local Green function of the correlated 
electrons is constructed, which is then used in the equation (|59|) to deliver the bath Green function Qa^abii^)- This 
matrix is the input to the impurity model. Solution of this model delivers the new self-energy Adabi^^) and the 
process is iterated towards self-consistency. 

Notice that once the DMFT self-consistency is reached, the process can either be stopped or continued since the 
Green function ^/^^(k, io;) delivers new charge density of the system which modifies the Hartree and LDA exchange- 
correlation potentials in the expression (|58|1 . In this respect, the LDA+DMFT method assumes a double iterational 
loop, the internal one over the self-energy and the external one over the density. This is precisely dictated by the 
spectral density functional stationarity condition. We illustrate such loop on Fig. 2| Note that in order to access 
accurate total energies and remove ambiguity that the LDA Green function (and not any other one) is used as an input 
to the DMFT calculation, this density self-consistency loop needs to be carried out. Our application to the volume 
expansion in Pu described later in this paper involves solution of the SDFT equations allowing the full relaxation of 
the charge density. 

Iterations over the density are not complicated to incorporate in the programs for electronic structure calculations. 
The spherical part of the density at a given site can be written approximately using the atomic sphere approximation 
as an integral over the partial density of states Ni{E) till the Fermi level Ep 

P{r) ^T. f ' Ni{E)vKr. E)dE (60) 

where (y9/(r, E) are the solutions of the radial Schrodinger equation. Usually these are constructed using spherical part 
of the LDA potential but in the present context the non-local self-energy operator needs to be utilized^: 



{-yli-E)^i{T,E)+ / Mefs{ry.E)^i{r',Eydr^O 



(61) 
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FIG. 4: Illustration of self-consistent cycle in spectral density functional theory with LDA+DMFT approximation: double 
iterational cycle consists of the innner DMFT loop and outer (density plus total energy) loop. 

Expression can be simplified further if we assume a Taylor expansion of (r, E)=ifi (r, i?^/ ) + {E — E^i)(pi {r, E,ji) 
around some linearization energies E^i taken at the centers of gravities of the occupied energy bands, i.e. Ei,i — 
jf^ENi{E)dE/ jf^Ni{E)dE. During the iterations, DMFT modifies the density of states, ANi{E), which 
leads to the change Ap(r) of the density. The latter has a feedback onto the change A(pi{r, E) [or changes 
Aipi{r, E^i), A(pi{r, E^i)]. If we assume that these changes are small comparing to the original LDA values, we 
can work out a linear response relationship for the density 

M^)-E ^Ni{E)^fir,E)dE + J2 Ni{E)^i{r,E)A^i{r,E)dE + Y,Ni{EFWf{r,Ep)AEp (62) 

and for the LDA potential 

AVLDAir) = / P^dr' + ^^Ap(r) (63) 
J \r-r'\ dp 

Thus, to first order, these are the quantities which are iterated in the external density loop of the spectral density 
functional as shown on Fig. 0] 

The main physical point of the LDA+DMFT methodology is identification of a subset of the correlated orbitals {Xa} 
which is separated from the full Hilbert space {x^}. In the case like Pu, this is the subset of f-electron orbitals. In 
other situations, this subset can be isolated based on physical grounds. If {Xa} is appropriately constructed Wannier 
representation, this subset may describe the bands crossing the Fermi level. We expect the dynamical self-energy 
corrections to appear at first place only within the subset {x^,}. However, changes in the electronic densities of states, 
ANi{E), will appear for all light and heavy electrons. 

We did not discuss so far the relaxation of the screened Coulomb interaction W(r, r', ilu), which, in principle, needs 
to be done during the self-consistency in parallel to Q{r, r', iuj). We stress that the short-range behavior appears only 
for the local effective susceptibility P{r, r', ioj) in and means its fast decay when r departs from r'. Contrary, the 
function W(r,r',iw) can be as long range as the bare Coulomb interaction if necessary. This is dictated by Eq. I|33() 
and is similar to the relationship H2U|I between Q and Mint- 

The locality assumption for V should simplify the self-consistency over W. This should be faster than the one 
employed in the full GW method which formally tries to compute full n(r,r',za'). In the language of local orbital 
representation {Xq}; this means computation of all matrix elements for nQ^(R — R',ia;) or its Fourier transform 
Ilai3(k,iuj) as compared to the site-diagonal {Srw) or small cluster cases of SDFT. This will be discussed below in 
connection to the recently proposecii2ii& local version of the GW method. 

So far we did not mention the problem of the optimal choice of the double counting corrections. This is discussed 
below in connection to the LDA+U method. 



18 



E. Double Counting and LDA+U Method 



Historically, the LDA+U method has been introduced^ as an extension of the local spin density approximation 
(LSDA) to treat the ordered phases of Mott insulating solids. In this respect it is a natural extension of LSDA. 
However, this method was first to recognize that a better energy functional can be constructed if not only the density 
but the density matrix of correlated orbitals is brought into the density functional. We have discussed the correlated 
subset {Xa} ^^'^ local Green functions Giocapi^^) in connection to the LDA+DMFT method. The density matrix 
Uab is related to the correlated subblock of the local Green function 

n,b = r^e'"°^G,o,,a6M (64) 

Therefore, the LDA+U method can be viewed as an approximation (Hartree-Fock approximation) to the spectral 
density functional within LDA+DMFT. 

The correct interaction energy among the correlated electrons can be written down explicitly using the Hartree-Fock 
approximation. In our language the LDA+DMFT interaction energy functional H55() is rewritten in the form 

^ LDA+U [p, n.ab] = Eh [p] + ^'^c'^ [p]+^u [nab] - ^DC [nc] (65) 
where the functional form <I>;7[riab] is known explicitly: 

^U[nab\ = ^ '^{Uacbd - Uacdb)nabncd (66) 
abed 

Here, indexes a, 5, c, d involve fixed angular momentum I of the heavy orbitals and run over magnetic m and spin a 
quantum numbers. The on-site Coulomb interaction matrix Uabcd is the on-site Coulomb interaction matrix element 
^a=ai3=b-^=c&=d appsared in H54() which is again taken for the subblock of the heavy orbitals. Note that sometimes 
Uabcd is defined as Vji^^^^^^^s^d- 

The double counting term '^Dc[nab] needs to be introduced since both the L(S)DA and U terms account for the 
same interaction energy between the correlated orbitals. This includes in first place the Hartree part. However, the 
precise form of the double counting is unclear due to non-linear nature of the LDA exchange-correlation energy. In 
practice, it was proposed- that the form for ^dc is 



^Dc[nc] = ^Un,{n, - 1) - -J[nl{nl - 1) + ni{ni - 1)]. (67) 

where U = ^^l+l)-^ J2ab Uabab , J ^ U ~ 21(21+1) Y^abi^abab " Uabba)a-nd where = J^aeh '^aaSa^a, flc = fil + 

Some other forms of the double countings have also been discussed in Ref. IHsl 

The minimization of the functional TLDA+u[p,nab] is now performed. The self-energy correction in appears 
as the orbital dependent correction Aiab — ^S>~^ '■ 

Mab = T = > (Uacbd - Uacdb)ncd (68) 

oriab ^ 

cd 

= ^ = ^abViUc -\)- J« - \)\ (69) 
OUab 2, I 

While the correction is static, it is best viewed as the Hartree-Fock approximation to the self-energy 7Wa(,(ia;) 
within the LDA+DMFT method. Note that such interpretation allows us to utilize double counting forms within 
LDA+DMFT as A^(r,r',zoo) or A^(r,r',zO). Note also that the solution of the impurity problem collapses in the 
LDA+U method since the self-energy is known analytically by formula Ht)8|) . 

From a practical point of view, despite the great success of the LDA+U theory in predicting materials properties 
of correlated solids^ there are obvious problems with this approach when applied to metals or to systems where the 
orbital symmetries are not broken. They stem from the well-known deficiencies of the Hartree-Fock approximation. 
The most noticeable is that it only describes spectra of magnetically ordered systems which have Hubbard bands. 
We have however argued that a correct treatment of the electronic structure of strongly correlated systems has to 
treat both Hubbard bands and quasiparticle bands on the same footing. Another problem occurs in the paramagnetic 
phase of Mott insulators: in the absence of any broken symmetry the LDA+U method reduces to the LDA, and the 
gap collapses. In systems like NiO where the gap is of the order of eV, but the Neel temperature is a few hundred 
Kelvin, it is unphysical to assume that the gap and the magnetic ordering are related. For this reason the LDA+U 
predicts magnetic order in cases that it is not observed, as, e.g., in the case of PupS. 
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F. Local GW Approximation 



We now discuss the relaxation of the screened Coulomb interaction W{r,r',iuj) which appeared in the spectral 
density functional formulation of the problem. Both LDA+DMFT and LDA+U methods parametrize the interaction 
W with optimally screened set of parameters, such, e.g, as the matrix Uabcd appeared in Ht)6|) . This matrix is supposed 
to be given by an external calculation such, e.g., as the constrained LDA method^^. To determine this interaction self 
consistently an additional self-consistency loop described by the equations and (|46|l has to be switched on together 
with calculation of the local susceptibility 7'(r,r',iw) by the impurity solver. This brings a truly self-consistent ab 
initio method without input parameters and the double counting problems. 

A simplified version of this method has been recently proposecii2ii& which is known as a local version of the GW 
method (LGW). Within the spectral density functional theory, this approximation appears as approximation to the 
functional sdf[Gioc, Wioc] taken in the form 

"^LGwiGioc, Wioc] = -^TrGiocWiocGioc (70) 

As a result, the susceptibility V{r,r',iuj) is approximated by the product of two local Green functions, i.e.T-" = 
—25"^ LGW 1 5Wioc — GiocGioc, and the exchange-correlation part of our mass operator is approximated by the local 
GW diagram, i.e. M^c = S'Slgw/SGioc = -GiocWioc- 

Thus, the impurity model is solved and the procedure can be made self-consistent: For a given Mint and V, the 
Dyson equations H20|l . H37|l for Q and W are solved. Then, the local quantities G/oc, Wioc are generated and used to 
find new Mint and V thus avoiding the computation of the bath Green function Qq after H45|) . and the interaction V, 
after g^. 

Note that since the local GW approximation (|70|l is relatively cheap from computational point of view, its imple- 
mentation on a cluster and for all orbitals should not be a problem. The results of the single-site approximation for 
the local quantities have been developed independently and reported in the literaturei^. The cluster extension is 
currently being performed and the results will be reported elsewhere^. 

Note finally that the local GW approximation is not the only one which can be implemented as the simplified 
impurity solver. For example, another popular approximation known as the fluctuational exchange approximation 
(FLEX) can be worked out along the same lines. Note also that the combination of the DMFT and full GW diagram 
has been recently proposed^^'^^ and a simplified implementation for Ni^-Z', and for a model Hamiltonian-^ have been 
carried out. This procedure incorporates full k-dependence of the self-energy known diagrammatically within GW 
together with the additional local DMFT diagrams. 



III. CALCULATION OF LOCAL GREEN FUNCTION 



The solution of the Dyson equations described in the previous section for a given strongly correlated material 
requires the calculation of the local Green function during the iterations towards self-consistency. This is very similar 
to the procedure in the density functional theory, when the charge density is computed. A big advantage of DFT is the 
use of Kohn-Sham orbitals which reduces the equation (|22|l for the Kohn-Sham Green function to a set of one-particle 
Schrodinger's like equations for the wave functions. As a result the kinetic energy contribution is calculated directly 
and the evaluation of the total energy of a solid is not a problem. Here, a similar algorithm will be described for the 
energy-dependent Dyson equation, the solution in terms of the linear-muffin-tin orbital basis set will be discussed, 
and the formula for evaluating the total energy will be given. 



A. Energy Resolved One-Particle Representation 



We introduced the auxiliary Green function C/(r,r',ia;) to deal with the kinetic part of the action in SDFT. It 
satisfies to the Dyson equation JSJ. Let us now introduce the representation of generalized energy-dependent one- 
particle states 

g{r,r ,iuj) = 2^—^ (71) 

g-\vy,tu) = ^^«^(r)(z^ + Ai-i?k,.)<-.(r') (72) 
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where the left V'kjw ^^"^ right tp^j^ (r) states satisfy to the following Dyson equations: 



-V^ + V;,t(r) + Vff (r)]^^.^(r) + J Al,,(r, r', zc^)7^£-^(r')cir' = i^kj^V^'-W (73) 
-V^ + Ve^tir') + yH(r')]<-^(r') + / ^^,^(r)X,e(r, r', zc^)rfr = E^.^^jt^Jr') (74) 



[we dropped the imaginary unit for simplicity in the notation V'kjw which shall be thought as a shortened version of 
Tpj^jir, ico)]. These equations should be considered as the eigenvalue problems with complex non-hcrmitian self-energy. 
As a result, the eigenvalues i^kjc^ [a shortened form for E]s^j{iLu)] being the same for both equations are complex in 
general. The explicit dependency on the frequency ilu in both eigenvectors and eigenvalues comes from the self-energy. 
Note that left and right eigenfunctions are orthonormal 

dr^ptjJr)i:^rJr)=S,,, (75) 

and can be used to evaluate the charge density of a given system using the Matsubara sum and the integral over the 
k-space: 

where 

We have cast the notation of spectral density theory in a form similar to DFT. The function gkjw is the Green function 
in the orthogonal left/right representation which plays a role of a "frequency dependent occupation number". 

It needs to be pointed out that the frequency dependent energy bands Ekju represent an auxiliary set of complex 
eigenvalues. These are not the true poles of the exact one-electron Green function G(r, r', z) considered at complex z 
plane. However, they are designed to reproduce the local spectral density of the system. Note also that these bands 
Ekjz are not the true poles of the auxiliary Green function t/(r,r',z). The latter ones still need to be located by 
solving a non-linear equation corresponding to the singularities in the expression (|71|l after analytic continuation to 
real frequency. For a one-band case this equation is simply: z + — — , whose solution delivers the quasiparticle 
dispersion Z^. General knowledge of the poles positions Z^j will allow us to write an alternative expression for Q 
which is similar to H71|l . but with the eigenvectors found at Z^j thus carrying out no auxiliary frequency dependence. 
These poles are the real one-electron excitational spectra in case ^ is a good approximation to G. However, the use 
of ifTT)) is advantageous, since it avoids additional search of poles and allows direct evaluation of the local spectral and 
charge densities the system. 

The energy-dependent representation allows us to obtain a very compact expression for the total energy. As we 
have argued, the entropy terms are more difficult to evaluate. However, they are generally small as long as we stay 
at low temperatures. The pure kinetic part of the free energy expressed via [see, Ea. (|39ll ] 

Tr\ng -TrMefjG ^Tj^e''^"^ J drdr' \ng{r,r' ,iuj) - tJ^ J drdr' Meffir,r' ,iuj)g{r' ,r,iuj) (78) 

iui icj 

needs to be separated onto the energy and entropy terms. Both contributions can be evaluated without a problem, 
but in light of neglecting the entropy correction in the interaction part, we concentrate on evaluating the kinetic 
energy only: 

r^e'-°" / dr [{-Vl)gir,r'M]^^^, ^ T^e'""" ^ ~ ^ (79) 

The SDFT total energy formula is now arrived by utilizing the relationship E\^^^ — (^kj^| — + -^e//|^kju;) = 
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EsDF = T^e*""^ ^5kj^£;kj^ - J drdr'Meff{r,r',tLj)g{r',r,iuj) + 

+ j dvV,^t{v)p{v) + ^JdrVH{r)p{r) + ^tJ^ J dvdv' M^^y ,t^)Gioc{v' ,v,iuj) (80) 

If the self-energy is considered as input to the iteration while the Green function is the output, near stationary point, 
it should have a convergency faster than the convergency in the Green function. 

It is instructive to consider the non-interactive limit when the self-energy represents a local energy-independent 
potential, say, the ground-state Kohn Sham potential of the density functional theory. This provides an important 
test for our many-body calculation. It is trivial to see that in the DFT limit, we obtain the Kohn- Sham eigenfunctions 



^ Vk,(r) (81) 
<-Jr) ^ Vk,(r) (82) 
-£'kjw — * -Ekj (83) 

and the one-electron energy bands are no longer frequency dependent. The sum over Matsubara frequencies in 
the expression for the charge density H76|l can be performed analytically using the expression for the Fermi-Diraq 
occupation numbers: 

and the formula (|76|) collapses to the standard expression for the density of non-interacting fermions. The total energy 
expression (|80|l is converted back to the DFT expression for the total energy since the eigenvalue E\^j^ becomes the 

DFT band structure E'kj, and the summation over Matsubara frequencies T X^i^j e^^^^gkjw gives according to H84() 
the Fermi-Diraq occupation number /^j . The standard DFT expression is recovered: 

Edft = /kj^kj - j drVeff{r)p{r) + J drK.t(r)p(r) + drV^^ (r)p(r) + E^^p] (85) 

where E^,, = (V'kjl - + Veffl^^j) = (V-kjl - + K.t + Vh + KdV^kj)- 



B. Use of Linear Muffln-Tin Orbitals 



The next problem is to solve the Dyson equation for the eigenvalues. The sophisticated basis sets developed to 
solve the one-electron Schrodinger equation can be directly used in this case. We utilize the LMTO method described 
extensively in the past literature^^iiiii as it provides a minimal atom-centered local orbital basis set ideally suited 
for the electronic structure calculation. Within the LMTO basis, the full Green function is represented as a sum 

g(r,r',z^) =^^X^(r)e„Mk,^^)xr(r') (86) 

k a/3 

and, as we have argued in the previous section, the matrix QapCk^iuj) needs to be considered as a variable in the 
spectral density functional. The stationarity yields the equation for the Green function 



Go^piKi^) = R«w + /i)0(k) -/i(")(k) -A^,„t(k,iw) 



-1 

a/3 



(87) 



where the overlap matrix OQ^(k) ~ (xa|x^),the non-interacting Hamiltonian matrix ^-^"^(k) — (Xq| — V'^-|-Vea;t(r)|x^) 
and the self-energy formally comes as a matrix element 

A1,„t,a/3(k,iw) = J drdr'xa*{Tc)Mint{r,r',iuj)xpir') (88) 
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over the LMTOs. Again, it is worth to point out that the self-energy here depends on k via the orbitals even if the 
single-impurity case is considered. In calculations performed on a cluster, the self-energy will also pick its non-trivial 
k-dependence coming from the nearest sites. 

While formally valid, the present approach is not very efhcient since the Green function Q{r,r',iuj) has to be 
evaluated via H86I) . This is the k-integral which has poles in a complex frequency plane, and integrating singular 
functions need to be performed with care. In this respect, we adopt the eigenvalue representation (|71|l . We expand 
the energy-dependent states in terms of the LMTO basis {x^} as follows: 



(89) 



The unknown coefficients Aq"''^'^^^'' are now the quantities which have to be considered as variables in the spectral 



density functional. The stationarity yields the equations 

{'k,iu!) - E]^j^Oaf3(k) 

E^a'"'"^ [/li°j(k) +X,;„t^„^(k,ZC^) - £;k,^OaMk) 



4^"^^^ = 



= 



(90) 
(91) 



These are the non-hermitian eigenvalue problems solved by standard numerical methods. The orthogonality condition 
involving the overlap matrix is 



(92) 



Q/3 



Note that the present algorithm just inverts the matrix H87|l with help of the "right" and "left" eigenstates. The 
Green function (|87|) in the basis of its eigenvectors becomes 



g^p (k, z^) = E IT- (93) 

This formula can be safely used to compute the Green function as the integral over the Brillouin zone, Eq. (jSBJ, 
because the energy denominator can be integrated analytically using the tetrahedron method^-'^ . 

Our next topic here is the evaluation of the bath Green function t/o(r, r', iw). It can be found from the integral 
equation 

go{r,r',iu;) = Gioc{r,r' ,iu;) - J dv" dv'" Gioc{vy ,iLo)M^nt{v" ,v"' ,iu)go{v"' ,v,iu) (94) 

where r and r' run over Q.ioc In order to solve this equation, it is useful to represent r = p + R, r' = p' + R', 
and redenote C/o(r,r',«w) = g^MR-'iP, p' ,'ii^),Gioc{r,r' ,iuj) = GiocMB' {p, p'^it^), Mint{r,r' ,iuj) = Mint,RR'{p-, p' 
Considering one atom per unit cell let us see how this can be solved using single-kappa LMTO-ASA method. The 
generalization to multiatomic systems with multiple kappa basis sets as well as inclusion of full potential terms in the 
calculation can be done along the same lines. The form of the LMTO basis function inside the sphere centered at R 
is 

X^(r) = X^(p)e^*'^ = e^'^^ic&f (p) + ^ *i (p)^L] (95) 

L 

where S*^^ are the structure constants of the LMTO method and where <I>^''^(p) are linear combinations of the 
solutions of the radial Schrodinger equation taken at spherical part of the potential and matched to the spherical 
Hankel (H) and Bessel (J) functions at the muffin-tin sphere boundary as well as their energy derivatives taken at 
some set of energies E^i at the center of interest. The local Green function can be represented in this basis set as 
follows 

Gi,,,nn'{p.p'M^Y. E *L^(P)GL"c':LL'fl/(*^)*i'^*(p') (96) 

LL' ^^u—H^J 
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where the matrices G^^^^^pji,{iu),Gf^l^j^^, j^,{iw),Gi"^^^^^ (indexes R and R' are restricted 

to a cluster) are given by the following Brillouin zone integrals 

Here are the original eigenvectors A^^'"-^'^) and A^^^^j^ = Ea^a^'^^^L' ^S^)"^ = Ea ^La^a'"^'^ are 

the convolutions of the eigenvectors with the LMTO structure constants. We now utilize a similar representation for 
the bath Green function 

where the matrices GQaRi3R'i'^'^)T^OaRL'R'i'''^)'So^R0R'i'^'^)'SoiRL'R'i'''^) *-an be found from the following Dyson 
equation (where the matrices sizes have been enlarged by a factor of 2) 

0L7hL'R' (*^) = G<i::l~RL'R' (*^) + -^i^r.LL'i^' (^^) (99) 

with the self-energy matrices are defined as follows 

■^£Ll'«'(*^) = / dpdp'<^>i^^*ip)M^r.tMR'ip,p\^u;)^^',\p') (100) 

The solution of the impurity model with ^o'iflL'i?' (*'^) delivers new matrix elements (|1UU|I . The k-dependent self- 
energy H88(l to be used in constructing the new Green function in H87|l is found first by restoring the k-dependence 

from the cluster ^^^^^'^^'(k, iuj) — J2r-R' ■^int^RL'R'i''''^)^^^^^'^ ^ and second, restoring the k-dependence of the 
LMTO basis as follows 

L L' LL' 

(101) 

In practical calculations performed with the LDA+DMFT method for Pu, only the subset of orbitals {Xa} treated 
as correlated (f electrons of Pu) and a single-impurity case is considered. It is useful to separate the Hartree and 
LDA exchange-correlations terms. Instead of dealing with the non-interacting Hamiltonian in (jH7J), we can rearrange 
the contributions to arrive 



GafiiK ii^) = + f^)0{k) - h^^^ik) ~ AM{k, iuj)] ^ (102) 

where h^J^^{k) = {x»\ + VLDA{r)\x^) with 14^^(0 = Kxt(r) + ^^(r) + F,^,^^(r). The matrix AA^„0(k, ic^) = 
Ma,0{k,iLu) - V^if{k), where 7W„^(k,iw) = S^aSpbMabiKitu) and V^fik) = S ^^aS pbV^^^ {k) represent the DMFT 
correction and double counting term described by Ht)9|l . These matrices are non-zero within the correlated subset. 
To accelerate the calculation of the impurity model, we can parametrize the self-energy matrix as Ai^J^^\iu)) — 
A^Ij^' (iaj)($i''' 1$^'^''). With such parametrization, the local Green function which enters the Dyson equation should be 

defined as follows G'ioc,af)(jt^) = Ekj E^i/ ^a(^)'^ ^'^^a^\^''h^)^(^)^ I {i^^iJ--E\^jS). This represents the generalization 
of a partial-density-of-state formula of the LMTO method. The bath Green function can be found from the equation: 
Cyo.ab(i^)^^ = Giocahij-^)^^ ^ -f^^ab ^^^) cau bc passed to the impurity solver. The latter delivers a new self-energy 
Ail^J (iuj) which is then multiplied by ($1^'' |$['^'') and used to reconstruct new k-dependent self-energy after (|l()l|l . 
Such procedure preserves all k-dependent information coming from the orbitals. 



IV. APPLICATIONS TO PLUTONIUM 



This section describes the application of the theory to Plutonium. Pu is known to be an anomalous metalS2.. It has 
six crystallographic structures. Starting from the low temperature a phase (0 to 100 C) with 16 atoms per unit cell 
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it shows a series of phase transitions and ends up in relatively simple fee 5 (300-450 C) and bee e phases (500-650 
C) just before it melts. The temperature dependence of atomic volume in Pu is anomalous^. It shows an enormous 
volume expansion between a and 6 phases which is about 25%. Within the 5 phase, the metal shows negative thermal 
expansion. Transition between S and higher-temperature e phase occurs with a 5% volume collapse. Also, Pu shows 
anomalous resistivity behavior^^ characteristic for the heavy fermion systems, but neither of its phases is magnetic. 
The susceptibility is small and relatively temperature independent. The photoemissionSi shows a strong narrow 
Kondo-like peak at the Fermi level consistent with large values of the linear specific heat coefhcient. 

Density functional based LDA and GGA calculations describe the properties of Pu incorrectly. They predict 
magnetic orderingS&. They underestimataSi equilibrium volume of the S and e phase by as much as 30% , the 
largest discrepancy known in LDA, which usually predicts the volume of solids within a few % accuracy even for 
such correlated systems as high temperature superconductors. Despite this, the volume of the a phase is predicted 
correctly by LDA^^. Since the transport and thermodynamic properties of a and S Pu are very similar, the nature of 
the a phase and the reason why LDA predicts accurately its structure and volume is by itself is another puzzle. 

To address these questions several approaches have been developed. The LDA+U method was applied to 5 PuffliSS. 
It is able to produce the correct volume of the S phase for values of the parameter U~A eV consistent with atomic 
spectral data and constrained density functional calculations. Similar calculation has been performed by a so-called 
orbitally ordered density functional method^. However, both methods predict Pu to be magnetic, which is not 
observed experimentally. The LDA+U method is unable to predict the correct excitation spectrum. Also, to recover 
the a phase within LDA+U the parameter U has to be set to zero which is inconsistent with its transport properties 
and with microscopic calculations of this parameter. Another approach proposed22i in the past is the constrained LDA 
approach in which some of the 5f electrons, are treated as core, while the remaining are allowed to participate in band 
formation. Results of the self-interaction-corrected LDA calculations have been reportedZi, as well as qualitative 
discussion of the bonding nature across the actinides series has been givenZ^. 

Thus, the problem of Pu is challenging because its f-electrons are close to the Mott transition^. It provides us a 
crucial test for our quantitative theory of strong correlations. A short version of this work has appeared alreadyS^. 
Our implementation is based on the self-consistent LDA+DMFT method and uses the LMTO method in its tight- 
binding representation^'^. Spin -orbit coupling effects arc important for actinide compounds and have been included 
in the calculation for Pu. The "full potential" terms have been neglected in the calculation through the use of the 
atomic sphere approximation with a one- kappa LMTO basis set. The necessary k-space integrals for evaluating Green 
functions and charge densities have been carried out using the tetrahedron method using (8,8,8) grid in the Brillouin 
zone. 

We study in detail the total energy as a function of the parameter U and give our predictions for the volumes in 
a, S, and s phases. A comparative analysis of the one-electron spectra in both a and S phases is also given and 
the comparison with the photoemission experiment^ is made. Since the dynamical mean field theory requires the 
solution of the Anderson impurity model for the multi-orbital f-shell of Pu, we have developed a method which, 
inspired by the success of the iterative perturbation theoryiS., interpolates the self-energy between small and large 
frequencies. At low frequencies, the exact value of the self-energy and its slope is extracted from the Friedel sum rule 
and from a slave-boson treatment^LS*^. This approach is accurate as it has been shown recently to give the exact 
critical value of U in the large degeneracy limit24 At high frequencies the self-energy behavior can be computed based 
on high-frequency moments expansionsMiiBi. The result of interpolation can be encoded into a simple form for the 
self-energy like the following continuous fraction expansion 

I](icj) = S(icx)) + — (103) 



where the unknown coefhcients are determined to satisfy known conditions in the low- and high-frequency limits. 
This kind of self-energy fits Quantum Monte Carlo (QMC) data in large region of parameters, such as U and doping, 
and where this comparison is at all possible (small degeneracy and high temperature). Thus, our approach interpolates 
between four major limits: small and large iw's valid for any U as well as small and large U's (band vs. atom) valid 
for any iuj. The analytical continuation to the real frequency axis is not a problem with the present method and avoids 
the use of the Pade^^ and maximum entropy'* based techniques. 

Complete details of the method can be found in Ref. |42. Here we only mention technicalities connected to the 
f-electrons of Pu where we deal with the impurity Green functions which are the matrices 14x14. However, for the 
relativistic f-level in cubic symmetry, the matrices can be reduced to 5x5 with 4 non-zero off-diagonal elements. The 
solution of such impurity problem is still a formidable numerical problem. We therefore make some simplifications. 
First, the off-diagonal elements are in general small and will be neglected. We are left with the 5f^/^ state split into 
2 levels which are 2-fold (Fy), and 4-fold (Fg) degenerate, and with the 5f''/^ state split into 3 levels which are 2-fold 
(Fg), 2-fold (F7) and 4-fold (Fg) degenerate. Second, since in Pu the intermultiplet spin-orbit splitting is much larger 
than the intramultiplet crystal field splitting (> 5:1), we reduce the problem of solving AIM for the levels separately 
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by treating the 5f^/^ and Fg levels as one 6-fold degenerate level, and the h'C 1"^ Fg, Fy and Fg levels as another 
8-fold degenerate level. 

A. Calculation of Volume 

We perform our calculations for 6 and e phases of Pu having simple foe and bcc structures respectively. Since our 
method does not yet allow us to treat complicated lattices, only a simplified study of the a phase which formally has 
16 atoms per unit-cell will be reported. The total energy as a function of volume is evaluated self-consistently using 
formula (|80() . It contains all the electrons including the core electrons. The local density approximation includes 
generalized gradient corrections after Ref. Since the LDA-I-DMFT approximation 1)55(1 is used, we subtract from 
the LDA the average interaction energy of the f electrons in the form H67|l of the double counting term and then add 
improved estimates of these quantities using the self-consistent solution of the impurity model. 

To illustrate the importance of correlations, we discuss the results for various strengths of the on-site Coulomb 
interaction U . Fig. [S] reports our theoretical predictions. First, the total energy as a function of volume of the 
fee lattice is computed. The temperature is fixed at 600-?^, i.e. in the vicinity of the region where the 8 phase is 
stable. C/=0 GGA curve indicates a minimum at 1//Vb=0.7. This volume is in fact close to the volume of the a 
phase. Certainly, we expect that correlations should be less important for the compressed lattice in general, but 
there is no sign whatever of the 8 phase in the \J—Q calculation. The total energy curve is dramatically different for 
[/>0. The details depend sensitively on the actual value of I] . The behavior at U —4 eV shows the possibility of a 
double minimum; it is actually realized for a slightly smaller value of U. We find that for U — 3.8 eV, the minimum 
occurs near V^/V5=0.80 which corresponds to the volume of the a phase if we allow for monoclinic distortions and 
a volume-dependent U. When U increases by 0.2 eV the minimum occurs at V/Vs = 1.05 which corresponds to the 
volume of the 6 phase, in close agreement with experiment. Since the energies are so similar, we may expect that as 
temperature decreases, the lattice undergoes a phase transition from the i5 phase to the a phase with the remarkable 
decrease of the volume by 25%. 
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FIG. 5: Total energy as a function of volume in Pu for different values of U calculated using the LDA-f DMFT approach. Data 
for the fee lattice are computed at T=600K, while data for the bcc lattice are given for T=900K. 



We repeated our calculations for the bcc structure using the temperature T — 900K where the e phase is stable. 
Fig. O shows these results for U = A eV with a location of the minimum at around V/Vs — l.Q3. While the theory 
has a residual inaccuracy in determining the 6 and e phase volumes by a few percent, a hint of volume decrease with 
the (5 — > e transition is clearly reproduced. Thus, our first-principles calculations reproduce the main features of the 
experimental phase diagram of Pu. 

Note that the values of U ~4 eV which are needed in our simulation to describe the a 6 transition, are in good 
agreement with the values of on-site Coulomb repulsion between f-electrons estimated by atomic spectral datcS, 
constrained density functional studies2&, and our previous LDA-I-U studies^. 
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The double-well behavior in the total energy curve is unprecedented in LDA or GGA based calculations but it is a 
natural consequence of the proximity to a Mott transition. Indeed, recent studies of model Hamiltonian systemsiSiS 
in the vicinity of the Mott transition show that two DMFT solutions which differ in their spectral distributions can 
coexist. It is very natural that allowing the density to relax in these conditions can give rise to the double minima as 
seen in Fig. |5l 

B. Calculation of Spectra 

We now report our calculated spectral density of states for the fee structure using the volume V/Vs=0.8 and 
V/Vs=1.05 corresponding to our theoretical studies of a and S phases. To compare the results of the dynamical 
mean-field calculations with the LDA method as well as with the experiment, we discuss the results presented in Fig. 
ini Fig. Ela) shows the density of states calculated using LDA+DMFT method in the vicinity of the Fermi level. Solid 
black line corresponds to the S phase and solid grey line corresponds to the a phase. We predict the appearance of 
a strong quasiparticle peak near the Fermi level which exists in the both phases. Also, the lower and upper Hubbard 
bands can be clearly distinguished in this plot. The width of the quasiparticle peak in the a phase is found to be 
larger by 30 per cent compared to the width in the 6 phase. This indicates that the low-temperature phase is more 
metallic, i.e. it has larger spectral weight in the quasiparticle peak and smaller weight in the Hubbard bands. Recent 
advances have allowed the experimental determination of these spectra, and our calculations are consistent with these 
measurements^^. Fig. IHI^b) shows the measured photoemission spectrum for 6 (black line) and a (gray line) Pu. We 
can clearly see a strong quasiparticle peak. Also a smaller peak located at 0.8 cV for the (5-phase can be found. We 
interpret it as the lower Hubbard band. 

The result of the local density approximation is shown on Fig. EJa) by dashed line. The LDA produces two peaks 
near the Fermi level corresponding to 5f^/^ and 5f^/^ states separated by the spin-orbit coupling. The Fermi level falls 
into the dip between these states and cannot reproduce the features seen in photoemission. We should also mention 
that LDA-I-U fails completelyiSSiSS, to reproduce the intensity of the /-states near the Fermi level as it pushes the 
/-band 2-3 eV below the Fermi energy. This is the picture expected from the static Hartree-Fock theory such as 
the LDA-I-U. Only full inclusion of the dynamic effects within the DMFT allows to account for both the quasiparticle 
resonance and the Hubbard satellites which explains all features of the photoemission spectrum in d Pu. 

The calculated by LDA-I-DMFT densities of states at Ef equal to 7 st./[eV*cell] are consistent with the measured 
values of the linear specific heat coefficient. We still find a residual discrepancy by about factor of 2 due to either 
inaccuracies of the present calculation or due to the electron-phonon interactions. However, these values represent an 
improvement as compared to the LDA calculations which appear to be 5 times smaller. Similar inaccuracy has been 
seen in the LDA-I-U calculation^. 

A simple physical explanation drawn from these studies suggests that in the S phase the f electrons are slightly on 
the localized side of the interaction-driven localization-delocalization transition with a sharp and narrow Kondo-like 
peak and well-defined upper and lower Hubbard bands. It therefore has the largest volume as has been found by 
previous LDA-I-U calculationsSMii which take into account Hubbard bands only. The low-temperature a phase is 
more metallic, i.e. it has larger spectral weight in the quasiparticle peak and smaller weight in the Hubbard bands. 
It will therefore have a much smaller volume that is eventually reproduced by LDA/GGA calculations which neglect 
both Coulomb renormalizations of quasiparticles and atomic multiplet structure. The delicate balance of the energies 
of the two minima may be the key to understanding the anomalous properties of Pu such as the great sensitivity 
to small amounts of impurities (which intuitively would raise the energy of the less symmetric monoclinic structure, 
thus stabilizing the 6 phase to lower temperature) and the negative thermal expansion. Notice however, that the a 
phase is not a weakly correlated phase: it is just slightly displaced towards the delocalized side of the localization- 
delocalization transition, relative to the S phase. This is a radical new viewpoint in the theoretical literature on Pu, 
which has traditionally regarded the a phase as well understood within LDA. However, the correlation viewpoint is 
consistent with a series of anomalous transport properties in the a phase reminiscent of heavy electron systems. For 
example, the resistivity of a-Pu around room temperature is anomalously large, temperature independent and above 
the Mott limiliSi (the maximum resistivity allowed to the conventional metal). Strong correlation anomalies are also 
evident in the thermoelectric power—. 

V. CONCLUSION 

In conclusion, this work describes a first principles method for calculating the electronic structure of materials where 
many-body correlation effects between the electrons are not small and cannot be neglected. It allows simultaneous 
evaluation of the total free energy and the local electronic spectral density. The approach is based on the effective 
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Energy relative Ep, eV 

FIG. 6: a) Comparison between calculated densities of states using the LDA+DMFT approach for fee Pu: the data for 
V jVi — 1.05, \J = 4.0 eV (black line) , the data for V jVi, = 0.80, V — 3.8 eV (gray line) which correspond to the volumes 
of the 5 and a phases respectively. The result of the GGA calculation (dotted line) at V jVi = 1 (C/ = 0) is also given, b) 
Measured photemission spectrum of & (black line) and a (grey line) Pu at the scale from -1.0 to 0.4 eV (after Ref. 16511 . 



action functional formulation of the free energy and is viewed as spectral density functional theory. An approximate 
form of the functional exploits a local dynamical mean-field theory of strongly correlated systems accurate in the 
situations when the self-energy is short ranged in a certain portion of space. The localization is defined with reference 
to some basis in Hilbert space. It does not necessarily imply localization in real space and is treated using a general 
basis set following the ideology of the cellular dynamical mean field theory. Further approximations of the theory, 
such as LDA+DMFT and local GW are discussed. Implementation of the method is described in terms of the 
energy-dependent one-particle states expanded via the linear mufhn-tin orbitals. Application of the method in its 
LDA+DMFT form is given to study the anomalous volume expansion in metallic Plutonium. We obtain equilibrium 
volume of the 8 phase in good agreement with experiment with no magnetic order imposed in the calculation. The 
calculated one-electron densities of states are consistent with the results of the photoemission. Our most recent 
studies^^ of the lattice dynamical properties of Pu address the problem of the 6 ^ e transitions and show good 
agreement with experimenljSi. 

Alternative developments of the LDA+DMFT approach by several groups around the world discuss other appli- 
cations of the dynamical mean field theory in electronic structure calculations. The results obtained are promising. 
Volume collapse transitions, materials near the Mott transition, systems with itinerant and local moments, as well as 
many other exciting problems are beginning to be explored using these methods. 
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